
****************************************************************
*paid sex:
  tab rel_sex, m 
  tab rel_sex, m nol  
  tab  REL_PAY_01 rel_sex, col ch
  tab  rel_sex_pay rel_sex, col ch

  preserve
drop if rel_sex==1
  tab  REL_PAY_01 sexe, col ch
    tab  REL_PAY_01 rel_sex_pay, r ch
    tab  rel_sex_pay sexe, col ch
restore


preserve
drop if PERDU_Q6_01==1
drop if rel_sex_Q6==1
tab REL_PAY_Q6_01 sexe, col ch
bysort sexe: tab REL_PAY_Q6_01 TEST_POS_ID, col ch
bysort sexe: tab REL_PAY_Q6_01 TEST_Q6_POS_ID, col ch
bysort sexe: tab REL_PAY_Q6_01 RISK_INFECTE_01, col ch 
tab  REL_PAY_Q6_01 rel_sex_pay_Q6, r ch
bysort sexe: ttest AGE, by(REL_PAY_Q6_01)
restore

preserve
drop if PERDU_Q12_01==1
drop if rel_sex_Q12==1
tab REL_PAY_Q12_01 sexe, col ch
restore

    tab con_profsex sexe, col ch
	tab CON_PROFSEX_01 sexe, col ch
	tab CON_PROFSEX_SPRE sexe, col ch
    tab con_profsex REL_PAY_01, col ch

preserve
drop if PERDU_Q6_01==1
tab CON_PROFSEX_01 sexe, col ch
restore

preserve
drop if PERDU_Q12_01==1
	tab CON_PROFSEX_01 sexe, col ch
restore

  *Q0
    preserve 
drop if rel_sex==1
  graph pie, over(REL_PAY_01) plabel(_all percent)
restore

  preserve 
drop if rel_sex==1
  graph pie, over(REL_PAY_01) by(sexe) plabel(_all percent)
restore

*Q6 : REL_PAY_Q6_01

tab rel_sex rel_sex_Q6, m
tab rel_sex_Q6, m
tab rel_sex_Q6, m nol

preserve
 drop if MONTH_Q6==.
 drop if rel_sex_Q6==1
   tab REL_PAY_Q6_01 sexe, col ch
   list sexe if REL_PAY_01==0 & REL_PAY_Q6_01==1
  restore

preserve
 drop if MONTH_Q6==.
 drop if rel_sex==1
   tab REL_PAY_Q6_01 sexe, col ch
  restore
  
    *veo entre los retenidos, cuÃ¡ntos habÃ­an dicho en basal paid sex
  preserve
 drop if MONTH_Q6==.
   tab REL_PAY_01 sexe, col ch
  restore
  
    preserve
 drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
   tab REL_PAY_01 sexe, col ch
  restore
  
   preserve 
drop if MONTH_Q6==.
 drop if rel_sex_Q6==1
  graph pie, over(REL_PAY_Q6_01) plabel(_all percent)
restore
  
       preserve 
drop if MONTH_Q6==.
 drop if rel_sex_Q6==1
  graph pie, over(REL_PAY_Q6_01) by(sexe) plabel(_all percent)
restore

*Q12
  preserve
 drop if MONTH_Q12==.
   tab REL_PAY_Q12_01 rel_sex, col ch
  restore

  preserve
 drop if MONTH_Q12==.
  tab rel_sex_Q12, m
    tab rel_sex_Q12, m nol
  restore

	preserve
 drop if MONTH_Q12==.
  drop if rel_sex_Q12==1
         tab  rel_sex_pay_Q12, m
   tab REL_PAY_Q12_01 sexe, col exact
       tab  REL_PAY_Q12_01 rel_sex_pay_Q12, r ch
     restore
  

   preserve 
drop if MONTH_Q12==.
  drop if rel_sex_Q12==1
  graph pie, over(REL_PAY_Q12_01) plabel(_all percent)
restore
  
   preserve 
drop if MONTH_Q12==.
  drop if rel_sex_Q12==1
  graph pie, over(REL_PAY_Q12_01) by(sexe) plabel(_all percent)
restore


****************logistic BIVARIANTE *******************************************************
* VA AUMENTANDO GRADUALMENTE PAID SEX CON LA EDAD
preserve 
drop if rel_sex==1
logistic REL_PAY_01 b2.sexe
logistic REL_PAY_01 b1.AGE_15_24
logistic REL_PAY_01 i.AGE_3cat
 restore 
 
tab AGE_15_24, m
tab AGE_15_24, m nol 


  *PARECE MÃS PAID SEX EN MENOS ESTUDIOS
*feb 2023. hay muy pocos con bajo nivel de estudios, pienso que no compensa meter en modelo
  tab ETUD_123, m
    tab ETUD_123, m nol
 tab REL_PAY_01 ETUD_N, col ch
  tab REL_PAY_01 ETUD_123, col ch
  tab REL_PAY_01 ETUD_3, col ch
    tab REL_PAY_01 ETUD_12, col ch
	
  tab PROF_3cat, m
 tab PROF ETUD_123, m


 preserve 
drop if rel_sex==1
     logistic REL_PAY_01 b2.ETUD_123
 restore
 
  preserve 
drop if rel_sex==1
       tab REL_PAY_01 PAR_ETUD, col ch
    tab REL_PAY_01 PAR_ETUD_23, col ch
 restore
 
      tab profess, m
	  tab profess, m nol
	  
	   preserve 
drop if profess==.   
      tab PROF_3cat, m
	  tab PROF_3cat, m nol
 restore
 
 preserve 
 drop if rel_sex==1
 drop if profess==.   
      logistic REL_PAY_01 i.PROF_3cat
  logistic REL_PAY_01 PROF_EST_TRAB
restore 
 *va aumentando ligeramente a mayor nivel econ
   tab ECON_123, m

 preserve 
drop if rel_sex==1    
tab REL_PAY_01 ECOLE_12, col ch
  logistic REL_PAY_01 i.ECON_123
    logistic REL_PAY_01 i.ECON_23
restore 



*MEDIA DIARIO
tab1 M_J_dia M_R_dia M_TV_dia M_int_dia M_phone_dia M_Rev_dia, m

tab1 M_J M_J_dia M_J_012, m
tab ECON_123 M_J_dia, col ch
tab ETUD_123 M_J_dia, col ch
tab ETUD_123 M_J_dia, col ch
tab prop_voit M_J_dia, col ch
tab prop_TV M_J_dia, col ch
tab prop_ord M_J_dia, col ch

tab1 M_J_dia M_J_012, m

 preserve 
drop if rel_sex==1
logistic REL_PAY_01 MEDIA_dia_3plus
logistic REL_PAY_01 MEDIA_dia_4plus
logistic REL_PAY_01 M_J_dia 
logistic REL_PAY_01 M_J_dia i.ECON_123
logistic REL_PAY_01 i.M_J_012
logistic REL_PAY_01 i.M_J_FREC
restore 
  
preserve 
drop if rel_sex==1
	    logistic REL_PAY_01 M_R_dia
    	    logistic REL_PAY_01 i.M_R_012
logistic REL_PAY_01 M_R_FREC
restore


  preserve 
drop if rel_sex==1
	logistic REL_PAY_01 M_TV_dia
	logistic REL_PAY_01 i.M_TV_012
		logistic REL_PAY_01 M_TV_FREC
    restore 
*internet parece que algo influye
  preserve 
drop if rel_sex==1
	logistic REL_PAY_01 M_int_dia   
	logistic REL_PAY_01 M_int_FREC
	logistic REL_PAY_01 i.M_int_012
restore 

tab M_int M_J, m

  preserve 
drop if rel_sex==1
	    logistic REL_PAY_01 M_phone_FREC
    logistic REL_PAY_01 M_phone_dia
		logistic REL_PAY_01 i.M_phone_012
restore 
*ORACION PROTEGE

tab RELIG_5cat_sinperd ECON_123, r ch

 preserve 
drop if rel_sex==1
tab REL_PAY_01 RELIG_3cat, m
logistic REL_PAY_01 i.RELIG_3cat
logistic REL_PAY_01 RELIG_CAT
logistic REL_PAY_01 RELIG_REVEIL
logistic REL_PAY_01 i.RELIG_5cat_sinperd
 tab REL_PAY_01 RELIG_CAT, col ch
   tab REL_PAY_01 RELIG_3cat, col ch
      tab REL_PAY_01 RELIG_3cat, col ch
	 tab REL_PAY_01 RELIG_5cat_sinperd, col ch
 restore 

 preserve 
drop if rel_sex==1
logistic REL_PAY_01 EGLISE_d
logistic REL_PAY_01 EGLISE_wk_d 
logistic REL_PAY_01 b1.EGLISE_wk_d 
logistic REL_PAY_01 PRIE_d
logistic REL_PAY_01 PRIE_wk_d 
logistic REL_PAY_01 b1.PRIE_wk_d 
  tab REL_PAY_01 relig_prie, col ch
  tab REL_PAY_01 PRIE_wk_d, col ch
    tab REL_PAY_01 PRIE_wk_d, col ch exact
    tab REL_PAY_01 PRIE_d, col ch
 tab REL_PAY_01 relig_eglise, col ch
  tab REL_PAY_01 EGLISE_wk_d, col ch
    tab REL_PAY_01 EGLISE_wk_d, col ch
 tab REL_PAY_01 EGLISE_d, col ch
restore 

 preserve 
drop if rel_sex==1
logistic REL_PAY_01 EGLISE_wk_d i.RELIG_3cat
logistic REL_PAY_01 PRIE_wk_d i.RELIG_3cat
restore 
   *CASADOS MÃS PAID SEX PERO EN REGRESIÃ“N PIERDE SIGNIF., ASÃ QUE ES MÃS EL EFECTO PAID SEX
   tab marie, m
   tab part_viv, m
      tab part_viv PART_VIV_01, m

	   preserve 
drop if rel_sex==1
drop if sexe==1
   logistic REL_PAY_01 ENCEINTE_01
   logistic REL_PAY_01 ENCINTE_FOIS4
restore 

   preserve 
drop if rel_sex==1
tab marie, m
tab maritale, m
    logistic REL_PAY_01 MARIE_01
    logistic REL_PAY_01 i.MARITALE_4cat
logistic REL_PAY_01 MARITAL_MONO_v2
	 logistic REL_PAY_01 PART_VIV_01
restore 


  ***ALCOHOL MUY ASOCIADO!
   preserve 
drop if rel_sex==1
  logistic REL_PAY_01 i.ALCOHOL_0123
  logistic REL_PAY_01 ALCOHOL_dia
  logistic REL_PAY_01 ALCOHOL_dia_sem
 logistic REL_PAY_01 ALCOHOL_SEX
  logistic REL_PAY_01 REL_FORCE_01
    logistic REL_PAY_01 REL_FORCE_PRACT_01
 restore
 
 ***ninguno condom siempre entre los de paid sex
 *importante sex under 15
 preserve 
drop if rel_sex==1
 logistic REL_PAY_01 SEX_AGE_U15
 logistic REL_PAY_01 SEX_AGE_U18
 restore
 
    preserve 
drop if rel_sex==1
 logistic REL_PAY_01 PART_M
 logistic REL_PAY_01 PART_M_6M_01
 logistic REL_PAY_01 PART_N_6M_MEME_01
 restore
 
  preserve 
drop if rel_sex==1
   logistic REL_PAY_01 i.CON_012
  logistic REL_PAY_01 CON_01
 logistic REL_PAY_01 CON_0
 logistic REL_PAY_01 CON_SPRE
 tab REL_PAY_01 CON_SPRE, m
  logistic REL_PAY_01 CON_PROFSEX_01
 restore
 
  preserve 
drop if rel_sex==1
    tab REL_PAY_01 CON_PROFSEX_01, r ch
	    tab REL_PAY_01 con_profsex, r ch
 restore
 
 
   preserve 
drop if rel_sex==1
  logistic REL_PAY_01 REL_ORAL_01
  logistic REL_PAY_01 REL_ANALE_01
  logistic REL_PAY_01 rel_sex_HMX
 restore
 

list AGE if sexe==2 & REL_PAY_01
sum AGE if sexe==2 & REL_PAY_01
 
 list AGE if sexe==1 & REL_PAY_01
sum AGE if sexe==1 & REL_PAY_01

list n_econ if sexe==2 & REL_PAY_01
list n_econ if sexe==1 & REL_PAY_01
list ETUD_123 if sexe==2 & REL_PAY_01
list ETUD_123 if sexe==1 & REL_PAY_01

list alcool if sexe==2 & REL_PAY_01
list alcool if sexe==1 & REL_PAY_01

tab ALCOHOL_5_01 if sexe==1 & REL_PAY_01

list ALCOHOL_SEX if sexe==2 & REL_PAY_01
list ALCOHOL_SEX if sexe==1 & REL_PAY_01

 tab TEST_POS_ID if sexe==1 & REL_PAY_01
 tab TEST_POS_ID if sexe==2 & REL_PAY_01

  tab TEST_POS if sexe==1 & REL_PAY_01
 tab TEST_POS if sexe==2 & REL_PAY_01
 
  tab test_hiv if sexe==1 & REL_PAY_01
 tab test_hiv if sexe==2 & REL_PAY_01
 
   tab RISK_INFECTE_01 if sexe==1 & REL_PAY_01
 tab RISK_INFECTE_01 if sexe==2 & REL_PAY_01
 
   tab con if sexe==2 & REL_PAY_01
  tab con if sexe==1 & REL_PAY_01

  *MISCONCEPTIONS
      preserve 
drop if rel_sex==1
 tab REL_PAY_01 SIDA_DIEU_SORC_SI, col ch
  tab REL_PAY_01 SIDA_SORC_SI, col ch
   logistic REL_PAY_01 SIDA_DIEU_SORC_SI
 restore
  
   *HIV
   preserve 
drop if rel_sex==1
logistic REL_PAY_01 i.test_hiv 
logistic REL_PAY_01 i.test_hiv sexe
logistic REL_PAY_01 TEST_POS
logistic REL_PAY_01 TEST_POS_ID
  restore
 
 tab TEST_AVANT, m
 sum TEST_FOIS
 tab1 sida_risque sida_infecte, m
    preserve 
drop if rel_sex==1
logistic REL_PAY_01 TEST_AVANT
logistic REL_PAY_01 TEST_FOIS
logistic REL_PAY_01 TEST_FOIS_2
logistic REL_PAY_01 TEST_FOIS_2plus
logistic REL_PAY_01 PAR_TEST_01
logistic REL_PAY_01 RISK_01
 logistic REL_PAY_01 i.RISK_012
 logistic REL_PAY_01 RISK_ALTO
logistic REL_PAY_01 RISK_INFECTE_01 
 restore
 
 tab ist_diag, m
   preserve 
drop if rel_sex==1
 logistic REL_PAY_01 IST_diag_01
   logistic REL_PAY_01 IST_ULCER_01
 restore
 
   preserve 
drop if rel_sex==1
   tab REL_PAY_01 sante, col ch
   tab REL_PAY_01 SANTE_45, col ch
      tab REL_PAY_01 SANTE_12, col ch
 restore
 

 
 *MODELO DEFINITIVO 09/03/2023:
 	*CON ANAL SEX EN VEZ DE ORAL SEX:
	   preserve 
drop if rel_sex==1 
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ANALE_01 TEST_POS_ID 
	bysort sexe: logistic REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ANALE_01 TEST_POS_ID 
	restore 
	
	preserve 
drop if rel_sex==1 
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01 TEST_POS_ID 
	bysort sexe: logistic REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01 TEST_POS_ID 
	restore

	
	   preserve 
drop if rel_sex==1 
	logistic TEST_POS_ID REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia EGLISE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	bysort sexe: logistic TEST_POS_ID REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia EGLISE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	restore
	
		   preserve 
drop if rel_sex==1 
	logistic TEST_POS REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia EGLISE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	bysort sexe: logistic TEST_POS REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia EGLISE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	restore
	
		   preserve 
drop if rel_sex==1 
	logistic TEST_POS REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	bysort sexe: logistic TEST_POS REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	restore
	
			   preserve 
drop if rel_sex==1 
	logistic IST_diag_01 REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	bysort sexe: logistic IST_diag_01 REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	restore
	
preserve 
drop if rel_sex==1 
	logistic REL_FORCE_01 REL_PAY_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	bysort sexe: logistic REL_FORCE_01 REL_PAY_01 b1.AGE_15_24 PROF_EST_TRAB M_J_dia PRIE_wk_d MARIE_01 ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ORAL_01  
	restore
	
	 *MODELO  16/02/2023:
   preserve 
drop if rel_sex==1 
    logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.PROF_3cat M_J_dia PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.PROF_3cat M_J_dia PRIE_wk_d PART_M_6M_01 SEX_AGE_U15 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ECON_23 M_J_dia PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.PROF_3cat M_J_dia PRIE_wk_d PART_N_6M_MEME_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ECON_23 M_J_dia PRIE_wk_d PART_N_6M_MEME_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	restore
	
*MODELO :
   preserve 
drop if rel_sex==1 
    logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_123 i.PROF_3cat M_J_dia PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
  logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ANALE_01
  logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 PROF_TRAB M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 i.PROF_3cat M_J_dia marie PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
   logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 i.PROF_3cat M_J_dia i.RELIG_5cat PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
       logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ECON_123 PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
	          logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ECON_123 PRIE_wk_d PART_N_6M_MEME_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
restore

   preserve 
drop if rel_sex==1 
    bysort sexe: logistic REL_PAY_01 b1.AGE_15_24 i.ETUD_123 i.PROF_3cat M_J_dia PRIE_wk_d PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 restore

   preserve 
drop if rel_sex==1 
  bysort AGE_15_24:logistic REL_PAY_01 b2.sexe i.ETUD_12 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
restore

 preserve 
drop if rel_sex==1 
 logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_12 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 logistic REL_PAY_01 b2.sexe age_n i.ETUD_12 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
logistic REL_PAY_01 b2.sexe b1.AGE_15_24 i.ETUD_123 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 logistic REL_PAY_01 b2.sexe AGE i.ETUD_123 i.PROF_3cat M_J_dia TEST_POS_ID PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
 logistic REL_PAY_01 b2.sexe b1.AGE_15_24  i.ETUD_123 i.PROF_3cat M_J_dia TEST_POS_ID PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
  logistic REL_PAY_01 b2.sexe b1.AGE_15_24 MARITAL_MONO_v2 EGLISE_wk_d i.ETUD_123 i.PROF_3cat M_J_dia M_int_FREC M_R_dia TEST_POS_ID PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
   logistic REL_PAY_01 b2.sexe i.AGE_DHS i.ETUD_123 i.PROF_3cat M_J_dia PART_M_6M_01 ALCOHOL_dia_sem REL_ORAL_01
   logistic REL_PAY_01 b2.sexe b1.AGE_15_24 ETUD_3 i.PROF_3cat M_J_dia PART_M_6M_01 CON_01 ALCOHOL_dia_sem REL_ORAL_01
restore


 preserve 
drop if MONTH_Q6==.
    bysort sexe: logistic REL_PAY_Q6_01 b1.AGE_15_24 M_J_dia PART_M_6M_01 ALCOHOL_dia_sem REL_ORAL_01
 restore

 
 ***ASOC CARACTS Q6---PAID SEX Q6
 tab1 prop_TV PROP_TV_01 prop_ord PROP_ORD_01 prop_voit PROP_VOIT_01 REL_PAY_Q6_01, m
 
 preserve 
drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
  logistic REL_PAY_Q6_01 PROP_TV_01
    logistic REL_PAY_Q6_01 PROP_ORD_01
  logistic REL_PAY_Q6_01 PROP_VOIT_01
  logistic REL_PAY_Q6_01  PRIE_Q6_wk_d
  logistic REL_PAY_Q6_01 REPAS_Q6_01
  logistic REL_PAY_Q6_01 PREN_PETIT_DIA_01
   logistic REL_PAY_Q6_01 ALCOHOL_Q6_d_wk
   logistic REL_PAY_Q6_01 RISK_INFECTE_Q6_01
   logistic REL_PAY_Q6_01 PART_M_Q6_01
     logistic REL_PAY_Q6_01 PART_M_6M_Q6_01
	 logistic REL_PAY_Q6_01 PART_N_6M_MEME_Q6_01
	 logistic REL_PAY_Q6_01 CON_Q6_01
	 logistic REL_PAY_Q6_01 REL_ORAL_Q6_01
	 logistic REL_PAY_Q6_01 REL_ANALE_Q6_01
 restore

 *asoc Q0---> pay q6
     preserve 
drop if MONTH_Q6==.
drop if rel_sex_Q6==1
  logistic REL_PAY_Q6_01 b2.sexe 
  logistic REL_PAY_Q6_01 b1.AGE_15_24
  logistic REL_PAY_Q6_01 b2.ETUD_123
  logistic REL_PAY_Q6_01 PROF_EST_TRAB
   logistic REL_PAY_Q6_01 MEDIA_dia_4plus
  logistic REL_PAY_Q6_01 M_J_dia
    logistic REL_PAY_Q6_01 M_int_FREC
	logistic REL_PAY_Q6_01 i.RELIG_3cat
  logistic REL_PAY_Q6_01 PRIE_wk_d
  logistic REL_PAY_Q6_01 MARIE_01
  logistic REL_PAY_Q6_01 ALCOHOL_dia_sem
  logistic REL_PAY_Q6_01 SEX_AGE_U15
    logistic REL_PAY_Q6_01 PART_M
  logistic REL_PAY_Q6_01 PART_M_6M_01
    logistic REL_PAY_Q6_01 PART_N_6M_MEME_01
	    logistic REL_PAY_Q6_01 CON_01
    logistic REL_PAY_Q6_01 REL_ORAL_01
	    logistic REL_PAY_Q6_01 REL_ANALE_01
  logistic REL_PAY_Q6_01 TEST_POS_ID
  logistic REL_PAY_Q6_01  TEST_AVANT_01
  logistic REL_PAY_Q6_01  RISK_INFECTE_01
    logistic REL_PAY_Q6_01  SIDA_DIEU_SORC_SI
   restore
	
    preserve 
drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
  logistic REL_PAY_Q6_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB i.RELIG_3cat ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ANALE_01  
	bysort sexe: logistic REL_PAY_Q6_01 b1.AGE_15_24 PROF_EST_TRAB i.RELIG_3cat ALCOHOL_dia_sem SEX_AGE_U15 PART_M_6M_01 REL_ANALE_01  
	restore
	
	    preserve 
drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
  logistic REL_PAY_Q6_01 b2.sexe b1.AGE_15_24 PROF_EST_TRAB ALCOHOL_Q6_d_wk PART_M_6M_Q6_01 REL_ORAL_Q6_01  RISK_INFECTE_Q6_01 TEST_POS_ID
	restore
	
	 preserve 
drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
  logistic REL_PAY_Q6_01 PROP_TV_01
    logistic REL_PAY_Q6_01 PROP_ORD_01
  logistic REL_PAY_Q6_01 PROP_VOIT_01
  logistic REL_PAY_Q6_01 REPAS_Q6_01
  logistic REL_PAY_Q6_01 PREN_PETIT_DIA_01
   logistic REL_PAY_Q6_01 PART_M_Q6_01
	 logistic REL_PAY_Q6_01 PART_N_6M_MEME_Q6_01
 restore
	
	
	
	 preserve 
drop if MONTH_Q6==.
  drop if rel_sex_Q6==1
  tab REL_PAY_Q6_01 sexe, m
     restore
	



